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Abstract 

The C. Neumann system describes a particle on the sphere S n under the influence 
of a potential that is a quadratic form. We study the case that the quadratic 
form has I + 1 distinct eigenvalues with multiplicity. Each group of m a equal 
eigenvalues gives rise to an 0(m CT )— symmetry in configuration space. The combined 
symmetry group G is a direct product of t + 1 such factors, and its cotangent 
lift has an Ad*-equivariant Momentum mapping. Regular reduction leads to the 
Rosochatius system on S , which has the same form as the Neumann system albeit 
for an additional effective potential. 

To understand how the reduced systems fit together we use singular reduction 
to construct an embedding of the reduced Poisson space T*S n /G into M 3 "" 3 . The 
global geometry is described, in particular the bundle structure that appears as a 
result of the superintegrability of the system. We show how the reduced Neumann 
system separates in elliptical-spherical co-ordinates. We derive the action variables 
and frequencies as complete hyperelliptic integrals of genus I. Finally we prove 
a convexity result for the image of the Casimir mapping restricted to the energy 
surface. 



1 Introduction 

The C. Neumann system is one of the few examples of an integrable system with n degrees 
of freedom, where n is an arbitrary positive integer. It describes a particle moving on the 
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sphere S n under the influence of a linear force. As C.G. Jacobi's student, Carl Neumann 
studied the case n = 2 in his thesis [27]. The general case has been beautifully described 
by Moser [23], [2H [25]. Knorrer JTB] elucidated the connection to the geodesic flow on 
ellipsoids, see also [31] . Later Knorrer [TTJ also showed that near hyperbolic critical 
values of the Integral mapping the Kolmogorov non- degeneracy condition holds so that 
KAM theory can be applied. Whether this is true for all regular values is still open. 
Devaney [7J noticed that the system has transversal homoclinic intersections while being 
Liouville integrable. Ratiu [28J realised that the system can be written as dynamics on 
the adjoint orbit of SO(n) x Sym(n) on the Lie algebra of this semi-direct product. The 
action variables and their Picard-Fuchs equations haven been derived in [8]. 

All these works assume the generic case that the quadratic potential ^(x; Ax), x 6 
R' t+1 of the system embedded in M. n+1 has n + 1 distinct eigenvalues a$ < a\ < . . . < 
a n . The present paper studies the degenerate case in which the eigenvalues a u are not 
all distinct. This modification leads to the Rosochatius system, an integrable (n + 1)- 
parameter deformation of the Neumann system, which we derive and analyse. The fact 
that the degenerate Neumann system is also integrable has been shown in |32j. This is 
non-trivial because the smooth integrals (fl6j), see below) of the non-degenerate system 
found by Uhlenbeck [23J become singular in this limit. So even though equal eigenvalues 
make the system simpler in the sense that global symmetry is introduced, a non-trivial 
transition takes place. 

The non-degenerate Neumann systems on S n with n > 1 does not admit globally 
defined continuous symmetries. The symmetry group of the generic Neumann system is 
the direct product Z^ 1 of n + 1 factors 0(1) = Z 2 and hence discrete. Correspondingly, 
there are no global smooth actions, i.e. integrals of motion that globally generate periodic 
flows. By the Liouville-Arnold theorem almost everywhere there exist local actions, but 
usually they cannot be extended to smooth global integrals. Hence locally there is a free 
T"-action (defined by the commuting flows of local actions Ij, cf. [B]), but globally there 
is not even a T 1 -action. Nevertheless, the system can be integrated by separation of 
variables. Explicit solutions can be derived in terms of ^-functions of genus n, see [21] . 

The degenerate Neumann system on S n admits a large global symmetry group, and 
hence possesses corresponding global actions. For each group of m a equal eigenvalues 
we have an 0(m a ) symmetry that acts by rotation of the corresponding group of co- 
ordinates and by cotangent lift on the momenta. In general, the potential has £+1 distinct 
eigenvalues bo < . . . < bi of multiplicities mo, . . . , m^. The corresponding symmetry group 
of the Neumann systems on S n is G = O(m ) x 0(mi) x . . . x 0(mg). In the generic 
case £ = n that all eigenvalues are distinct we recover G = Z 2 +1 . The symmetry group G 
describes the freedom in the choice of the (non-unique) co-ordinate system in which A is 
diagonal. 

In most of the paper we consider the orthogonal groups as the relevant factors in the 
symmetry group G, because this allows a uniform treatment of the factors 0(m) no matter 
what m is. However, sometimes it is more appropriate to restrict in one or several factors 
to the special orthogonal group. This is the case when the singularities of the reduced 
phase space can be avoided, and it occurs where a factor of G is 0(2). The reason for this 
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is a subtle difference between the invariants of the 0(2) and the SO (2) actions, which is 
not present when the factor is 0(m) with m > 2. This difference becomes particularly 
important when studying Hamiltonian monodromy, which appears e.g. when there are 
two eigenvalues of multiplicity two each. The corresponding situation in the geodesic flow 
on the ellipsoid has been studied in [51 [6j. 

Our main result is the reduction of the degenerate Neumann system with the symmetry 
group G and a detailed description of the relation of the reduced system with £ degrees 
of freedom with the Rosochatius system. After describing the setting in the next section, 
and a review of orthogonal group actions in Section [21 this reduction is performed in 
Section HI When the momentum is weakly regular, Marsden-Weinstein reduction [20] 
leads to a system on S with an additional effective potential. Since not all momentum 
values are weakly regular we perform singular reduction, which for singular values leads 
to reduced phase spaces with singularities and simultaneously embeds all reduced phase 
spaces coming from (weakly) regular values as symplectic leaves of a Poisson structure in 
R 3e+3 . When m a > 3 for at least one index a the system is superintegrable, with a non- 
commutative symmetry group, and the lower dimensional tori are generically parametrised 
by a product of oriented Grafimannian manifolds Gr m<7 ,2- Furthermore we compute the 
relative equilibria and characterise the Energy-Casimir mapping. In Section we show 
how the Rosochatius system obtained by reduction of the degenerate Neumann system 
can be separated in elliptical-spherical co-ordinates. From this the non-trivial actions 
and frequencies are derived. In the final Section [6] we show that the set of critical values 
of the Integral mapping of the reduced system (consisting of the I independent integrals 
of motion that are not obtained from the symmetry group G) is topologically a 2^-tant 
when all Casimirs are nonzero. Finally, we prove that the image of the Casimir mapping 
restricted to the energy surface is a convex set, which in the limit of large energy tends 
to a convex polyhedron. 



2 The Neumann System 

The Neumann system describes a particle moving on a sphere under the influence of a 
quadratic potential. We use x = (x , X\, . . . , x n ) as co-ordinates in M n+1 in which the 
sphere S n is embedded as C\ = (x; x) = 1 with the standard Euclidean scalar product. 
The kinetic energy is T(x) = ~(x;x) and the potential is V(x) = |(x;Ax). Newton's 
equations are 

x = - W + Ax (1) 

where the Lagrange multiplier A gives the strength A = 2V — 2T of the constraining 
normal force. The initial conditions must be chosen such that C\ = 2(x;x) = 0. To 
preserve C\ we furthermore require C\ = 0, which determines A. By a rotation of the 
co-ordinate system we can always achieve that A is diagonal with diagonal entries a u , 
v = 0,1, ... ,n. The classical Neumann system has n = 2, and distinct eigenvalues of A 
ordered as clq < a\ < a 2 . We call a Neumann system (with arbitrary n) degenerate if A 
has multiple eigenvalues. 
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2.1 Global Hamiltonian description 

The canonical momenta are y = (jjq, yi, ■ ■ ■ y n )- They satisfy the canonical Poisson bracket 
relations [x u ,y K ] = 5 UK , [x u ,x K ] = [y v , y K ] = 0. Besides C\ — 1 the momentum vectors 
are constrained to be tangent to the sphere, C2 = (x; y) = 0. This embeds the phase 
space T*S n into IR 2n+2 . Since we are dealing with a constrained system for which x are not 
generalized co-ordinates the Lagrangian is degenerate and the Legendre transformation to 
the Hamiltonian does not work. Instead we just write down the Hamiltonian in canonical 
variables of the embedding space 

ff(x >y ) = i<y;y) + y(x). (2) 

To ensure that the constraints C\ — 1 and C2 = are respected by the corresponding 
Hamiltonian vector field, we modify the canonical Poisson bracket [..,..] to the Dirac 
Poisson bracket {..,..}, cf. e.g. [4]. We then confirm directly that {..,..} together with (T5]) 
lead to Newton's equations as well, see §5§ below. 

For the Dirac bracket one first calculates the bracket of the constraints. The only 
nonzero bracket in our case is [CijCy = 2C±. Alltogether they form the matrix 7^- = 
[Ci,Cj]. The original bracket is modified by [/, Ci][Cj, <7](7 _1 )y- This gives 

{/,<?} = If, 9) + ^- i [f,C 1 ][C 2 ,g] - JLfcCftCug] . (3) 

Sometimes we need this bracket in spaces of different dimensions, in which case we indicate 
it by a subscript, whence the above bracket reads {.. , ..}2n+2- This new Poisson structure 
{..,..} has C\ and Ci as Casimirs and is explicitly given by 

{x u ,x K } = 0, {x v ,y K } = 5 VK - , {y v> y K } = ■ ( 4 ) 

Here L UK = x u y K — x R y v are the components of the angular momentum J : y*]R n + 1 — > 
so(n + 1)* ~ R™(«+ 1 )/ 2 . The Poisson bracket {L UK ,L^} is nonzero if and only if exactly 
one of the indices u, k coincides with one of the indices /x, A. Using L VK = —L KU the single 
coinciding index can be moved so that we obtain all nonzero bracket relations from 

{L UK , L K \] = L\ u . 

Note that {L UK ,L K \} = [L UK ,L K \] because [Lij,Ci] = 0. Now Hamilton's equations 
/ = {/, H} restricted to T*S n are 

x = y, y = -VU+((x;VV)-2T)x. (5) 

Here we used the Casimir C 2 = and the identity ^2,L vn y K = x u YIVk- These equations 
are equivalent to Newton's equations ([T]). 
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The Neumann system with distinct eigenvalues a v is Liouville integrable. The poly- 
nomial integrals F , . . . , F n given in [23] are 

n j 2 

F v (x,y) + (6) 

and they are independent up to the relation = C\. When restricting to T*S n the 

Hamiltonian can be expressed as 

1 ™ 

Surprisingly, these integrals were not classically known [23] . 

Equations (jSJ) are valid for dynamics on S n with arbitrary potential. Integrable systems 
in this class of Hamiltonians include not only the Neumann system with V^(x) = ~ a v x t 
but also the spherical pendulum with V(x) = x or the Rosochatius system [221 EH EES] 
with V(x) = \ a u x 2 +w v /x 2 . For n = 2 and the linear potential we obtain the classical 
spherical pendulum on S 2 , see e.g. [3]. For n = 2 with the quadratic potential we obtain 
the classical Neumann system on S 2 [27]; the degenerate case with a\ = a-i ("quadratic 
spherical pendulum" ) has been studied in [31 E] • 

In this paper we are concerned with the cases in which A has multiple eigenvalues, so 
that there is additional rotational symmetry among the groups of co-ordinate axes x v with 
the same a v . Let b a , a = 0, . . . , £, denote the values of different coefficients a u where b a has 
multiplicity m„. We arrange groups of equal a v to be labeled consecutively and assume 
they are sorted by size bo < . . . < bp. Generally we use greek indices v, k, . . . ranging from 
to n and a, r, . . . ranging from to £, while latin indices i,k, . . . run through the index 
sets I a that contain the m a indices which have the same coefficient b a = for i £ I a . 
Then the potential can be written as 

V = b V + b 1 V 1 + ... + b e Ve, V a = -J2*l 

The potential is invariant under the symmetry group 

G = O(m ) x 0{mi) x . . . x O(m^) 

where 0(m a ) acts on the co-ordinates summed over in V a : the sum of their squares 
remains constant. The group G describes the ambiguity for the choice of an orthonormal 
basis of eigenvectors of A. Indeed, the co-ordinate system chosen above is obviously not 
unique, but any G-action on it gives another admissible co-ordinate system (possibly 
with opposite orientation) in which the Hamiltonian assumes the same form. Under 
the rotations from G the velocities transform (by tangent lift) in the same way as the 
co-ordinates, so the kinetic energy remains constant as well. 
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The integrals ([6]) become singular in the degenerate case. The integrals in the degen- 
erate case were found in [32]. For each set of m a equal coefficients aj = b a , % 6 I a the 
integrals are the angular momenta where i,k G I a and the function 




^ = £*. 2 +£££^V 



(7) 



Assembling the angular momenta into 



£ Ll 



we can express the Hamiltonian on T*S n as 



H 



^Y^b a F a + W a . 



Considering the case with equal coefficients the limit (denoted by a — > b 

below) of the generic case with all distinct the integrals F u of the generic case diverge. 
However, the following combinations remain finite, 



2.2 Basic Dynamics 

Some simple observations can be made directly from the differential equations of the 
Neumann system. Understanding this local information helps to put together the global 
picture later on. 

Observe that the differential equations ([T]) have the simple form x = Cx with the 
diagonal matrix C = diag((A — a v )). This may appear to be linear, but in fact A = 2V — 2T 
depends on all the variables. Still, for any choice of subset of d indices / = {vi, . . . , 
the 2d dimensional invariant subsystem x v = y u = 0, v ^ I, is again a Neumann system 
of lower dimension. 

In the simplest case d = 1 the invariant subsystem consists of two equilibrium points at 
(x, y) = (±e^, 0). The Hessian of the potential V in local co-ordinates (which are x with 
x v omitted) near the equilibrium has eigenvalues 2{a K — a u ), k/v. Denote the number of 
a K > a u by m" and the number of a K < a u by m s v . Then the spectrum of the equilibrium 
at ±e u has m" nonzero real eigenvalues and ml nonzero imaginary eigenvalues. For 
the coefficients b a with multiplicity m a there is a sphere S m " y ~ 1 of equilibrium points of 
dimension m a — l given by V a = |. All of them are equilibrium points of ([!]) because A — b a 





keia 



and the relation Ylt=o F a — C\ 



1 on S n is inherited. 



6 



exactly balances the force for all nonzero co-ordinates. These non-isolated equilibrium 
points are reflected by the number m a — 1 of zero eigenvalues. 

For d = 2 there is a Neumann subsystem with 1 degree of freedom in the plane with 
the chosen indices. The potential on the circle obtained by intersecting the plane with 
the sphere C% = 1 is periodic, and either has two minima and two maxima (like "twice" 
a pendulum), or is constant when the coefficients a v in this plane are equal. The types 
of motion inside this plane are easy to obtain, but the normal stability of these orbits is 
already non-trivial. These n(n + l)/2 families of periodic orbits in particular contain the 
non-linear normal modes of the n + 1 equilibria. 

For d = 3 we find a Neumann subsystem with 2 degrees of freedom. The generic 
motion takes place on 2-tori T 2 when the coefficients of the three planes are not all equal. 
In this case there exist isolated periodic orbits on the energy shell which are not inside an 
invariant plane with d = 2. When all three coefficients are equal the potential is constant 
on the sphere, and the system defines the geodesic flow on S 2 . In that case all the orbits 
are closed, and in fact great circles on S 2 . 

The hierarchy of invariant subsystems continues up to d — n. Each invariant set has 
a normal stability, and invariant manifolds associated to that normal stability. Together 
with the invariant tori they form an intricate system of invariant submanifolds that to- 
gether foliate the phase space. The foliation for the non-degenerate Neumann system has 
been studied in [8]. Interestingly, the foliation becomes more complicated when additional 
symmetry is introduced. This additional structure is described in this paper. 

Now we focus on those invariant subspaces for which the above index set I is equal to 
an index set l a of equal coefficients, given by 

{(x,y)6TT:x, = |/„ = 0V^} . (8) 

Consider an initial condition (x , y ) with y ^ inside this invariant subspace. Choose 
a new orthogonal co-ordinate system which has among its unit vectors the directions of 
xo and yo- This is always possible by an orthogonal transformation because xo and yo 
are orthogonal. Moreover, the 2nd order differential equation (pQ) is invariant under such 
a rotation, because we only need to rotate within the subset of co-ordinates l a for which 
there are equal coefficients. We already noted that the symmetry group G originates 
from this freedom of choice of co-ordinates. After this co-ordinate transformation, only 
two of the rricr co-ordinates are involved in the dynamics, while the remaining m a — 2 
co-ordinates remain constant. In this sense locally nothing new happens when m CT > 3. 
However, globally, the geometry does become more interesting for m a > 3: a different 
initial condition gives a solution that differs by a rotation. The set of all such solutions 
is given by the set of all oriented 2-planes in M m,T . This is a GraBmannian manifold and 
it appears as a coadjoint orbit below. 

Interestingly, most of the above analysis can still be applied even when the initial 
condition (x , yo) is not inside an invariant subspace (JBj). We then have co-ordinates 
in different l a and all these can again be simplified by a rotation, so that only two 
co-ordinates in each subsystem are non-constant. The difference is that the Lagrange 
multiplier in ([!]) depends on all the co-ordinates and momenta, and therefore is not a 
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constant anymore. In this way, after the co-ordinate transformation all the remaining 
equations are coupled in a non-trivial way. 

When at least one m a > 3 the system becomes superintegrable, cf . [2^1 E21 EH EEU [12] , 
having more integrals than degrees of freedom. This makes the phase space a ramified 
torus bundle, with regular fibres T e+e+1 . Here I + 1 is the number of distinct eigenvalues 
of A and £+1 is the number of groups of multiple eigenvalues. If m a < 2 for all a = 0, . . . , £ 
then £ + £ + 1 = n and the regular fibres of the ramified torus bundle are Lagrangian 
tori. Under regular reduction the fibres T e+e+1 are mapped to — the generic reduced 
motion takes place on £-tori. Hence £ is the number of frequencies needed to describe a 
generic reduced motion. For each of the £+1 groups I a with m a > 2 the full system turns 
out to have one additional frequency. 

In the following we are going to make the above more precise and describe the asso- 
ciated global geometry. Our treatment of the degenerate Neumann system heavily relies 
on the reduction with respect to the symmetry group 

G = O(m ) x ... x 0{mt) 

with subgroup 

G = O(m ) x ... x 0(m £ ~) , G/G £ l£ l 

whence we now recall some facts about the cotangent lifted group action of the orthogonal 
group 0(m). We shall later see that understanding this group action is the key to an 
efficient description of the symmetry reduction in the degenerate Neumann system. The 
following section also serves as a quick introduction to (singular) reduction and to the 
geometry of superintegrable systems. 

3 Orthogonal Group Action 

Consider the orthogonal group 0(m) of invertible m x m matrices g which preserve the 
Euclidean scalar product in R m , so that for x G R m we have (gx; gx) = (x;x). Let the 
group G = 0(m) act on T*R m = R 2m with co-ordinates (x, y) and symplectic structure 
Q = dx A dy. The rotation g e 0(m) gives a point transformation by left multiplication, 
the induced action on the cotangent bundel T*R m is obtained by cotangent lift. Since 
{q 1 ) 1 = 9 we fi n d the symplectic action 



Hence, the orthogonal group acts by rotating co-ordinates and conjugate momenta in the 
same way. The general ^-invariant Hamiltonian is a function of the invariants of this 
group action, whence it can be expressed as a function 



$ : 0(m) x T*W 
(&x,y) 



m 




(9) 



tf(x,y) = ff«x;x),(y;y),<x;y)) 



(10) 
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of the three basic invariants of the action <3>. In other words, these three basic invariants 
form a Hilbert basis, and since there are no constraining syzygies between these invariants 
this is a free Hilbert basis. A standard example is the motion in a central force field, for 
which the Hamiltonian reads 

tf(x,y) = i(y;y} + [/«x;x» 

with radial potential U. 

3.1 Momentum mapping 

By Noether's theorem each one-parameter symmetry leads to a conserved quantity. For 
multi-parameter groups like 0(m) the integrals of motion can be assembled into the 
momentum mapping of the group action. 

Lemma 3.1. Form > 2 the action of G = 0(m) on T*IR m given by $ 5 (x, y) = (gx,gy) 
has an Ad* -equivariant Momentum mapping given by 

J : T*R m — ► so(m)* 
(x, y) i-)- x <S> y — y <g> x . 

Any ^/i G so(m)* is a weakly regular value, which is regular for m = 2, 3. 

We recall that Marsden-Weinstein reduction [20J works for weakly regular momentum 
values. A value of a different iable mapping is called weakly regular if the tangent space 
at every point of the pre-image is given by the kernel of the Jacobian at that point. 
Lemma 13.11 is true in general for faithful cotangent lifted matrix Lie group actions, see, 
e.g., 12.1.1 in [21]. For the convenience of the reader and in order to introduce our notation 
we give a direct proof. 

Proof. Let £ G so(m) be the infinitesimal generator of g = exp£. The symmetry $ g is 
the time one mapping of the flow of the vector field = (£x, £y). This vector field has 
the Hamiltonian = (£x; y). The Momentum mapping J : T*W m — > so(m)* is defined 
by 

<J(x,y),0« = ^( x >y) , 

where the pairing (A,B) S0 = |trAB* is used in so(m). Using tr(xy') = (x; y) equa- 
tion (fTTj) is easily verified. The components of J are of course the ordinary angular 
momenta = Xiyt — Xkyi, and generates a rotation in the i/c-plane. 

When the symmetry $ 9 is applied to (x, y) the corresponding image of the Momentum 
mapping is transformed into 

(J o $ s )(x, y) = {gx)(ffyy - (ffy)(yx)* = s(xy* - yxV . (12) 
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This in turn is just the (co)adjoint action of the group on the (dual of the) algebra and 
reads Ad*-i/x = gfig" 1 for the orthogonal group. Hence J is Ad*-equivariant, 

Ad*_i oJ = Jo$ r 

The value fj, — is not weakly regular because J _1 (0) is not a smooth manifold (instead, 
it is a cone). For /i ^ the set J _1 (/i) is a submanifold of T*W m , and the tanget space 

7(x ,y„)J^) = {(«x + /3y 0) 7x - «y ) G T (x()iyo) (T*M m ) :a,ft 7 6l} 

is identical to the kernel of -DJ(x ,y ) at any point (x ,y ) G J~ 1 ( / u). When m = 2, 3 
the rank of DJ is full when /i^O, and so the value is regular. □ 

3.2 Group orbits 

All symmetry transformations applied to a point p G T*R m produce the group orbit 
G ■ p = {$ 9 (p) : g G G} of p. It is isomorphic to the quotient of the full group G by the 
isotropy subgroup G p = {g G G : & g p = p}; the 'motions' that fix p have to be divided 
out. For convenience we define O(0) to be one point. 

Lemma 3.2. Let p = (x, y) G J _1 (//) ; £/ien £/ie G-orbit of p under the action $ is 
isomorphic to 

G/Gp = 0(m)/0(m — d) = V m ^ (Stiefel manifold) 
where d = dimspanjx, y}. 

The Stiefel manifold V m ^ = 0(m)/0(m — k) is the set of /c-dimensional orthonormal 
frames in M. m . 

Proof. The isotropy groups G p are easily determined. There are three cases, 

• d — 2 : x Ify =^ G p = 0(m — 2): the generic case, 

• d — 1 : x || y =^> G p = 0(m — 1): the singular case, 

• <i = 0:x = y = ^> G p = 0(m): the trivial case. 

In the generic case the 'motion' takes place in the 2-dimensional plane spanned by x 
and y. The isotropy group is the subgroup of rotations that fixes this plane. In the 
singular case the 'motion' takes place on a line, and G p is the subgroup of rotations that 
fix this line. The group orbits of points p in these three cases are diffeomorphic to 

• d = 2 : G/G p = 0(m)/0(m - 2) = V m , 2 : Stiefel manifold, diml^ = 2m - 3, 

• d = 1 : G/Gp = 0(m)/0(m - 1) = V m>1 = S m - 1 : Sphere, dimS™" 1 = m - 1, 

• d = : G/Gp = {id}: Point (x, y) = (0, 0). 
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□ 



Now fix the momentum p = 3{p). When p' = & g (p) runs through the orbit G ■ p of p then 
p' = J(p') = Ad*-i(/x) by equivariance. The set of all such momenta is the group orbit 
G ■ p of p under the coadjoint action Ad*. It is isomorphic to G/G^ where the isotropy 
subgroup is given by G M = {g G G : Ad*_i/i = p}. 

Lemma 3.3. The G-orbit of p ^ under the coadjoint action Ad*p = gpg~ l is isomor- 
phic to 

G/G^ = 0(m)/(SO(2) x 0(m - 2)) = G m , 2 (Grafimann manifold). 

The Grafimann manifold G m ^ = 0(m)/(SO(2) x 0(m — k)) is the set of fc-dimensional 
oriented subspaces of IR' m and can alternatively be defined as SO{m) / (50(2) x SO(m—k)). 

Proof. First consider the case m = 2. Let e = ( ^ and for a G 0(2) require 

ae = ea. This implies that actually a G 50(2). For m > 3 note that any nonzero rank 2 
antisymmetric matrix can be brought into the form diag(e, 0). Hence it is enough to 
consider this special p = diag(e, 0) G so(m)*. The invariance gpg 1 = p then requires, 
writing g in block form, 

fa b\ fe 0\ _ fe 0\ fa b\ 
\c d) \0 OJ ~ \0 OJ \c d) ' 

In addition to a G 50(2) this gives b = c = 0. Since g G O(m) the other block d G 
0(m — 2). Conversely, any g of the form diag(50(2), 0(m — 2)) leaves p = diag(e, 0) 
invariant, hence G M = 50(2) x 0(m — 2). The coadjoint orbit G ■ p = G/G^ is the 
Grafimannian G mj2 = 0(m) / (50(2) x 0(m — 2)), the set of oriented 2-planes of R m . □ 

Let if be a ^-invariant Hamiltonian, H o $ 9 = H for any g G 0(m). Then by Noether's 
theorem the flow y?^ of the Hamiltonian leaves J invariant, J o ip l H = J. Hence, each 
component of J is a constant of motion. All the components = XiUk — XkUi are 
integrals of motion, but only 2m — 3 of them are independent, rank J = 2m — 3. Even 
when all components are independent (as for m — 3) they do not all commute. 

As an example consider the general Hamiltonian ffTUl) with m = 3, or the particular 
case of the motion in a central force field. As usual for m = 3 the dual of the Lie algebra 
so(3)* is identified with M. 3 . The generic initial condition (x, y) satisfies x x y ^ so 
that the motion takes place in the plane through the origin perpendicular to the angular 
momentum x x y. This plane is fixed by the reflection g G G p = 0(1). The group orbit 
of a point p = (x, y) is 0(3)/G p = 50(3). The angular momentum p is not constant 
along this group orbit. It is only constant when g G G^ = 50(2) x 0(1), i.e. when we 
fix the axis given by the angular momentum. When other rotations act on p the angular 
momentum takes on all possible values in the coadjoint orbit 0(3)/G M = 5 2 . This is the 
space of oriented planes in R 3 . In the singular case the angular momentum vanishes and 
the line of motion is fixed by rotations about this line, while in the trivial case the origin 
is fixed by every rotation. In both cases the coadjoint orbit is not the correct notion for 
the analysis of the geometry since p = 0, see the singular reduction below. 
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3.3 Regular reduction 

The regular reduction procedure is to fix a weakly regular momentum value /i6§* and 
factor out the G^-action from the set J -1 (/z) of points that have fixed momentum \i and 
form the quotient 

Pp = J- 1 (jj.)/G ll . 

For an abelian group we have = G, while in the non-commutative case only dynamics 
associated to group elements g whose coadjoint action preserves \i is factored out, hence 
G^. The reduced symplectic manifold P M then is the space of G^-orbits in J _1 (/i). The 
dimension of the reduced space is given by 

dim P^ = dim P - 2 dim G/G p + dim G/G^ 

where P is the original phase space. Thus, for G = 0(m) and P = T*W m we find 

dimP M = 2m - 2(2m - 3) + (2m - 4) = 2 . 

The regular reduction for the 0(m)-action $ is made effective by giving concrete models 
for the abstract reduced phase spaces. 

Lemma 3.4. The reduced space P M for fi ^ is diffeomorphic to the open half plane in 
M. 2 with global symplectic co-ordinates 



f=||x|[>0, 77=(x;y) 
so that = 1. The ^-invariant Hamiltonian (flUI) becomes 

,P 

e 

where the total angular momentum J = ||J|| * = VW is given by 
W = UrJJ t = J2 L i 



ff((,ri 1 J) = C^ + - 1 ^ ) (13) 



2 

Ki<k<m 



2 

J ik 



Proof. Since /i ^ we have £ > 0. The half plane £ > is of course diffeomorphic (even 
symplectomorphic) to M 2 . The co-ordinates (£,?/) are symplectic since 



^ d£ drj 1 ^ 2 



i=i 



The identity 

W = ||x|| 2 ||y|| 2 -(x;y) 2 = (||y|| 2 -r / 2 )£ 2 (14) 

allows to eliminate ||y|| 2 from H and drop it to H on P M . The fibre of the reduction 
mapping from the submanifold of constant total angular momentum j, 

J- 1 (j) = {(x,y)€T*R m : J(x,y)=j} 

to (£, 77) is the Stiefel manifold found before. □ 
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For m = 2 the above reduction is equivalent to introducing the radius of polar co-ordinates 
in the plane. For fixed nonzero angular momentum radius zero is impossible, the line £ = 
is not part of the reduced space, and the reduction is regular. 

Remark 3.5. The function W is a polynomial in the basic invariants, as opposed to 
J = y/W. The square root is not differentiable at and this difference between J and W 
becomes important in the singular reduction later on. From an algebraic point of view W 
is more fundamental, but from the symplectic point of view J is more fundamental, since 
J defines an action outside /i = 0, i.e. J is a Hamiltonian with 27r-periodic flow when 
/i ^ 0. In case m = 2 one may choose to reduce only the S'0(2)-action instead of the 
full 0(2)-action. This yields a fourth basic invariant L 12 and the Hilbert basis ceases to 
be free as the basic invariants become related by the syzygy L\ 2 + (x;y) 2 = | |x| | 2 | |y| | 2 , 
in particular W = L\ 2 . The reason m = 2 is special is that J can be redefined as 
J = J = L12, which is polynomial and provides a smooth passage through the critical 
value [A — 0. 

Since we have many more integrals (namely m(m— 1)/2, of which 2m — 3 are independent) 
the Hamiltonian H is superintegrable for m > 3. For a weakly regular value of the Energy- 
Momentum mapping 

(H, J) : T*R m — >Rx so(m)* (15) 

the pre- image is a two-dimensional set. It is a two-torus when compact, and there are two 
actions associated to it. One action is the total angular momentum J, the other comes 
from the reduced Hamiltonian. All the other 2m — 4 dimensions do not take part in the 
dynamics. The semi- local situation is described by the following result. 

Lemma 3.6. (Generalized action-angle co-ordinates) The symplectic form on T*M. m lo- 
cally near a weakly regular value of the Energy-Momentum mapping (fT5j) decomposes 
into 

n = fi M + dip a d j + n G .^ 

where f2 M = d£ A dt] is the symplectic structure on the reduced space P^, if) is conjugate to 
the action J, and Qg-^ is the symplectic structure on the coadjoint orbit space G m ^- The 
reduced Hamiltonian ( !T3|) only depends on the co-ordinates (£, 77, J). For a generic invari- 
ant Hamiltonian a system of generalized action-angle co-ordinates in the sense of I2h^ is 
obtained by the action-angle co-ordinates of the reduced one- degree- of- freedom system in 
(£,77) and by (tp, J) together with any set of symplectic co-ordinates on G m ^- 

Proof. Let us first of all check that J is an action, i.e. that its flow is 27r-periodic. Note 
that using (fH|) 

\v y J 2 = | |x| | 2 y - (x; y)x and ^V X J 2 = | |y | | 2 x - (x; y)y . 
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For each component Hamilton's equations for J therefore are 

ft) —ft)' M = K-S'y> Nl ?> )" 

Note that 5M is an involution, (S'M) 2 = —id, which follows from det M = 1, see (Il4p . 
Hence the eigenvalues of SM are ±i, and the 27r-periodic flow is 

= (id cos t + sin t) 

for every z = 1, ... ,m. The projection of the orbit to configuration space (or momentum 
space) is a circle in the plane spanned by the initial x and y, with radius ||x|| (or ||y||). 
This follows because (— (x;y)x+ ||x|| 2 y)/J is orthogonal to x and with length ||x||, 
similarly for the second row of SM and y. The flow of J yields the angle ip conjugate 
to J. The computation of the 2nd action depends on the Hamiltonian. 

Since the reduced Hamiltonian has one degree of freedom and depends on the action 
J all the other co-ordinates belong to the symplectic manifold of co-ordinates that do 
not appear at all in the Hamiltonian. This symplectic manifold of ignorable symplectic 
co-ordinates is given by removing the periodic flow of J from the fibre of the reduction 
mapping, Vm^/S 1 = G rn ^- Hence, the coadjoint orbit G ■ /i = G m ^, with its symplectic 
structure (see, e.g., |21j) is this symplectic manifold. □ 

Generalized action-angle co-ordinates in the sense of [22] are obtained from this construc- 
tion when if is a generic 0(m)-invariant function, so that the generic motion has two 
frequencies. For special H almost every orbit is periodic, and in this case the construction 
needs to be modified. 

In the picure developed in [TTJ [T2] there is a two-dimensional meadow of actions (which 
globally can be replaced by the image of the Energy-Casimir mapping, see below), on 
which there are flowers whose petals are two-tori parametrised by the angles conjugate to 
the actions, and the centre of the flower carrying the petals is the coadjoint orbit G ■ fi. 
More precisely, the pre-image 

^ = J-\j)= U J -1 ^)' j = IHI 

fi€Gfi 

is mapped onto the centre of the flower, J (J-}) = G m ^- Therefore the flower is a bundle 
over G mj 2, and following Nekhoroshev [22] the fibre of this mapping is T 2 . Since the 
Nekhoroshev actions are only locally defined, some flowers are still missing in the global 
picture. 

3.4 Singular reduction 

All we have said so far is only valid for weakly regular values. To understand what 
happens near \x = we need to use singular reduction. Note that setting J = in 
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the reduced Hamiltonian H ([TBI does correctly describe reduced motion on a line in 
configuration space. However, if fails to describe how the neighbouring reduced systems 
with two degrees of freedom and \i ^ limit to this special case. 

Lemma 3.7. The reduced phase space P^ is the image of the momentum level set J -1 (/i) 
under the Hilbert mapping 

X : T*R m — > R 3 
(x,y) H- (y,T,S) 

where we put V = |||x|| 2 , T = ^||y|| 2 and S = (x;y). When \i ^ then P^ ~ R 2 ; given 
by the sheet V > 0, T > of the two-sheeted hyperboloid 

2VT - l -S 2 = \j\ (16) 

where j = ||/i|| * For — it is the half cone 2VT = \S 2 , with V > 0, T > 0. The fibres 
of x ore the Stiefel manifolds V m ^ with d = 2 in the regular case and d = 1 or d = 
for fj, — 0. The reduced phase space P^, \i ^ fixed, is the symplectic leaf of the Poisson 
structure on M 3 with co-ordinates (V,T,S) satisfying the bracket relations 

[V, T] = S, [V, S] = 2V, [T, S] = -2T 

o/sI(2,R). The function W = J 2 on T*R m is a Casimir of this bracket. The reduced 
Hamiltonian H is a function ofV, T and S only, reading H = T + U(V). 

Proof. For the computation of the brackets it is enough to compute e.g., [x 2 , Xiyj\ = 2x 2 , 
etc. The induced bracket on R 3 with co-ordinates (S, T, V) (for which we use the same 
notation [..,..]) inherits the Jacobi identity. While one may straightforwardly compute as 
well that W is a Casimir with respect to this bracket, this also follows from W : T*R m — > 
R factoring through 50 (m)* by means of the Ad*-invariant expression W = L 2 k . The 
symplectic leaves are given by the quadratic form fll6p . with the given relations. For 
W — the reduced space is not a smooth manifold, but half of a cone (in general such a 
reduced phase space is a semi-algebraic variety). □ 

Remark 3.8. The singular reduction for \x = leads to a reduced system on a non- 
smooth manifold, with a singular point at the tip of the cone 4VT = 5* 2 . An alternative 
description of this singular reduced phase space is obtained by first restricting to any two- 
dimensional subspace that is invariant under the dynamics of H, e.g. the (xi, 7/i)-plane 
defined by xi = yi = for i > 2. The O(m) -act ion $ has a residual Z 2 -action on this 
plane given by $(7r) : (xx,yi) (— X\, —yi). Therefore the singular reduced phase space 
can also be viewed as R 2 /Z 2 . This is e.g. the half-plane y\ > with the boundary yi — 
identified with itself by means of (xi,0) ~ (— xi,0), which again gives a half cone. The 
fixed point of the residual Z 2 -action generated by $(7r) is the origin x\ — y\ = and 
projects to the singular point of the reduced phase space. 
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The relation to the regular reduction in Lemma [3.41 is that (£,77) = (y2V, S/\^2V) are 
symplectic co-ordinates on any regular symplectic leaf. It is, however, not a co-ordinate 
system for W — 0. When fi = then W = and x || y for every p G J _1 (0). Hence 
the fibre over the cone Pq is V m \ = S ,T ™~ 1 . There is no dynamics on this sphere. When 
m = 2 this S l is part of a regular invariant two-dimensional torus for almost all values 
of the energy. For m > 3 this is not possible since in this case no if?" 1-1 — bundle over S 1 
is homeomorphic to an m-torus. The upshot is that monodromy can be defined when 
m = 2, but not when m > 3, see below for more details. 

Combining the previous results we can now describe the global meadow of actions, 
together with its maximal dynamical tori and symplectic manifolds in which no dynamics 
takes place in the regular and singular cases. 

Lemma 3.9. (Energy-Casimir Mapping) The Energy-Casimir mapping 

EC = {H, W) : T*R m — ► R 2 

classifies the dynamics. When it has full rank the fibre is a T 2 -bundle over G m ^ (or a 
T 1 -bundle over V m p)- Where the rank is 1 we have a relative equilibrium and the fibre 
over this point is V m ,2, which is a T 1 -bundle over G m ,2- In both cases the periodic flow 
of J gives an S 1 -action such that the centre of the flower is Vm^/S 1 = G m ^- For \i = 
the fibre of EC~ x {h^ 0) is itself a ramified sphere bundle over the base space 

I (V, T, S) e R 3 : H(V,T,S) = h,2VT = \S 2 , V > 0,T > o} 

with regular fibres V m> i = S*" 1-1 and singular fibre V m fl (a single point) over (V,T, S) = 0. 

Proof. This follows from Lemmas 13.21 and 13.31 □ 

Note that we cannot say that the fibre over a regular point is T 2 x G m> 2, since the S 1 - 
bundle V m ^ does not posses a global section. In fact, for m = 3 it is the Hopf fibration of 
V 3 , 2 = SO(3). 

One way to view this classification result is to consider the two stage process of first 
considering the Energy-Momentum mapping and then the mapping of sums of squares 
from momentum space to R that defines the dynamically relevant Casimir. 

In T*M? all this comes down to a familiar picture. The situation is simplified by 
identifying so(3)* with M 3 . The orbit G • fi = G/G^ = G 3 ,2 = is the sphere with 
radius j = ||/x|| in momentum space. Already most of the group is involved in generating 
the coadjoint orbit, while only the S 1 of rotations about the axis \l actually generates 
dynamics. The motion in the configuration space M. 3 takes place in the plane orthogonal 
to fi, which is the projection of the set J -1 (/i). The flow of J generates an S 1 in this 
plane, and factoring this out of J _1 (yu) gives the reduced space R 2 . The set of all such 
planes is the coadjoint orbit. 

For the critical value G so (3)* = R 3 the sphere shrinks to the point G ■ 0. The group 
orbit G ■ p C T*R 3 turns into the sphere V^i = S 2 encoding the common direction of 
x || y. For (x, y) = this sphere shrinks to a point as well. 
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In higher dimension the regular situation is similar, but more complicated by the fact 
that there are additional relations between the angular momenta that define the coadjoint 
orbit. These syzygies, together with J = const define the coadjoint orbit G ■ \i = G mt2 
as a submanifold of momentum space. E.g. for m = 4 the single additional syzygy is the 
Plucker relation Li 2 L 3 4 + L14L23 + L13L24 = 0. Such identities hold for any 4 indices when 
m > 5, but all such identities with terms LijLki, i < j,k < I are not independent of each 
other on their common zero level set. As before G mj 2 is the space of all 2-dimensional 
oriented planes in M m . No dynamics takes place in this space. When J > the dynamics 
in the invariant plane is described by (£, rj) and the angle if) conjugate to J. 

4 Reduction of the degenerate Neumann System 

The results of the last section are now applied to the action of a direct product of orthogo- 
nal groups. Many statements just go through because they hold separately for each factor 
of the group. The essential difference is that we are not starting from the symplectic man- 
ifold T*M n with the standard symplectic structure [..,..], but from the Poisson manifold 
T*R n+1 with the Poisson structure {..,..} that has T*S n among its symplectic leaves. 
The most degenerate case of the Neumann system appears when all spring constants in 
the potential are equal. This is the case treated first. 

4.1 0(n+ 1) symmetry on S n 

When all coefficients in the Neumann system are equal the potential is constant because 
it is proportional to the Casimir C\ of the Poisson structure 01]). Hence, the Hamiltonian 
describes the geodesic flow on S n with the induced metric from M n+1 . 

Surprisingly little changes as compared to section |3] when systems on S n embedded in 
M n+1 are considered. The reason is that the constraining Casimirs are both invariants of 
the group action, in fact they are given by C\ = 2V and C 2 = S. 

Actually things are much nicer for this action on T*S n , since it is "almost" free. The 
reason is that the origin x = is no longer part of the configuration space. Moreover, 
since x 1 y on T*S n q y*]g>n+i only way to achieve x || y is to have y = 0, i.e. no 
momentum at all. If we consider sufficiently high energy (or geodesic flows with nonzero 
energy from the start) this is impossible. The following theorem is well known, since it 
describes the geodesic flow on the sphere; we have merely formulated it in the general 
framework developed here for the degenerate Neumann system. 

Theorem 4.1. (Dirac Bracket and S n ) Endowing T*M. n+1 with the Dirac bracket 
the additional Casimirs V = \C\ and S = C 2 are among the three basic invariants V, 
T , S and the reduced brackets vanish identically. The reduced space merely is the point 
H = T = J 2 /2. The only dynamics in phase space is that generated by J, it is motion 
along the great circle defined by the intersection of S n and the plane spanned by x and 
y. The quotient of the energy surface {(x, y) G T*S n : H = h > 0} by the flow of H is 
J-\V2h)/S l = G n+ia . 
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Proof. Brackets involving Casimirs are always zero, and {T, T} = as well. When H is 
fixed then J is also fixed and if H > then J assumes regular values. So the meadow for 
the flowers is just a single point when we fix the energy. For a geodesic flow a change of the 
energy merely re-parametrises the orbits. The petals of the flower are the one- dimensional 
tori parametrised by the angle conjugate to J. The symplectic leaves are given by the 
quotient of the energy surface by the action of J. For H = 1/2 the energy surface is the 
unit sphere bundle 

T*S n = {(x,y) G T*R n+1 : ||x|| = 1, (x; y) = 0, ||y|| = 1} . 

The flow generated by J with respect to the Dirac bracket is a rotation in the plane 
spanned by x and y, as before (the original symplectic bracket [J, Cj] with the Casimirs 
vanishes). Hence, the orbit is an (oriented) great circle. The set of all oriented great 
circles is the Grafimann manifold G n+1|2 , the set of oriented 2-planes in whence 

T-* cm / q1 /~i 
l O / O — (Jn+1,2 • 

The reversing symmetry (x, y) h-> (x, — y) gives another orbit. These are two orbits 
corresponding to each 2-plane, one for each orientation. □ 

The flow of H is the flow of J, up to scaling time by a constant. This is why the above 
result can also be obtained from Lemma 13.31 



4.2 The general degenerate case 

Our main result is concerned with the more complicated case in which neither all coeffi- 
cients of the potential are the same, nor all of them are different. For each group of m a 
equal coefficients there is an 0{m a ) symmetry group acting on the space spanned by the 
axes of the equal coefficients. In fact the action is an 0(m CT )-action on M m<T , as described 
in Section El Since the axes of the groups of equal coefficients are different, the group 
action of the direct product of the 0{m a ) is simply given by the direct product of their 
actions. The difference to Section [3] (besides the fact that now we have a direct product of 
orthogonal groups) is that we need to consider the Dirac bracket {..,..} on M n+1 instead 
of the standard symplectic structure [..,..]. 

The reduction with respect to the action of the joint symmetry group is regular when 
the momentum of each component is nonzero. The total angular momentum for each 
group of equal coefficients is an Ad*-invariant function. Such functions become Casimirs 
after dividing out the group. In addition there is a relation between these Casimirs that 
comes from the Casimirs C\ and Ci resulting from the original embedding of the sphere in 
M n+1 . Hence, we get a reduced system that has i degrees of freedom: one for each group 
of equal coefficients minus one for the constraint to be on the sphere. The reduction is 
regular when all the fixed momenta are nonzero. In this case the reduction leads to a 
simple system with an effective potential. When some or all of the momenta are zero 
singular reduction needs to be used. In this case the reduced system is embedded in a 
higher dimensional Euclidean space, and the reduced phase space is no longer a smooth 
manifold, though still a semi-algebraic variety. 
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Proposition 4.2. (Symmetry Group Action) The Neumann system with £+1 groups of 
equal coefficients with multiplicity m a > 1, a — 0, . . . , £ is invariant under the symmetry 
group 

G = O(m ) x 0(mi) x ■ ■ ■ x 0{me) 
with action $ on T*IR n+1 , ^ m a = n + 1, 
$ : G x T*R n+l — > T*R n+1 

(x,y)) (gx,gy)- 

Here g is a block diagonal matrix from G C 0(n + 1). 

Proof. This is clear from the structure of the degenerate Neumann system. Each group 
of m a equal coefficients b a = a i: i G I a admits an 0{m a ) symmetry acting on the co- 
ordinates Xi with i G l a and the corresponding momenta The sets l a are disjoint, 
hence the action is a direct product. □ 

When the potential of the Neumann system is written in general form as | (x; Ax) we see 
that the symmetry group is in fact the group that describes the ambiguity of the choice 
of a diagonalizing co-ordinate system when A has multiple eigenvalues. 

Theorem 4.3. (Momentum Mapping) Denote by G the subgroup of G containing the fac- 
tors 0(m a ) with m a > 2 and denote the number of such factors by £ + 1. The Momentum 
mapping J of the $ -action of G C G is the direct product of the Momentum mappings 
of the individual factors for which m a > 2. It is Ad* -equivariant with momentum /iG§. 
The isotropy groups are the direct products of the isotropy groups of Lemma \3.1\ To each 
partial mapping corresponding to indices I a with m a > 2 there is a Casimir W a given by 
the Ad* -invariant function 

i<k£l a 

Proof. Since each partial action of 0(m a ) is as described in Lemma 13.11 this result im- 
mediately follows. When m a = 1 the corresponding group only has a discrete factor 
0(1), which acts by reflection in configuration space leading to the 7r-rotation (x a ,y a ) i-> 
(— x a , —y a )- This discrete symmetry does not appear in the Momentum mapping. Hence, 
the momentum fi only contains the components for each group with m a > 2. □ 

Remark 4.4. When m a = 2 it is also useful to consider SO (2) reduction instead of the 
0(2) reduction. The reason is that then L^, i < k G I a is among the invariants and 
W a = L 2 ik can be replaced by J a = L^. In this way the action J a becomes a signed 
quantity that is differentiable also at = 0, because no square root needs to be 

extracted. This is particularly important at \\/J* a \\ = because then the fixed set of the 
discrete symmetry is accessible for the dynamics. It is not possible to define such signed J a 
when m a > 3 because then the plane of rotation cannot be defined when = 0, while 
for m a = 2 there only is one plane. 
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Denote by V„ the projection onto the co-ordinates with indices I a , so that "Po-x gives all 
the co-ordinates Xi with i e I a , similarly V a y gives the corresponding momenta. In each 
such subspace we have the standard Eulidean scalar product and norm. 

Theorem 4.5. (Regular Reduction) Let each component be nonzero, hence the 

Casimir W a ^ for each m a > 2. For each such a define reduced co-ordinates by 

\\P<T~X\\ 

When rricr = 1 define £ CT = x a and T] a = y G . The reduction mapping R from (x, y) to (£, 77) 
is a Poisson mapping, i.e. it preserves the Dirac bracket {..,..}. The G-reduced phase 
space is an open subset ofT*S e . The reduced Hamiltonian reads 

H = li2vl + v„ ^ = ij>£ + |F- ( 17 ) 

CT=0 (T=0 ^< 7 

where formally W a = when m a = 1. 

Proof. The verification that this is a Poisson mapping can be done by direct computation. 
The reduction is as described in Lemma 13.41 now applied separately to each group of 
indices I a with m a > 2. The fact that the reduction does not lead to all of S e comes from 
the fact that £ CT = is impossible when \ \fJL a \\ 7^ 0. Nevertheless, the global Casimirs C\ = 
^2 — 1 and C 2 = Yl £.o-Va = restrict the motion to the reduced phase space T*S e . □ 



Clearly, the reduced Hamiltonian f Tl7|) defines an integrable system on S . This system 
was first studied by Rosochatius who separated it in elliptical-spherical co-ordinates. 

Proposition 4.6. (Generalized action-angle co-ordinates) The Poisson structure (j4j) on 
T*M n+1 locally near a weakly regular value of the Energy-Momentum mapping (H, J) has 
the bracket relations 



{£*,rir} = S aT -^, a,r = 0,... 

{Va,Vr} = kfcgk*, a,r = 0,... 



{L ik ,L kl } = L H , i, k,l e I a , a = 0, . . . ,£ 

with all further bracket relations being zero. The reduced Hamiltonian (TIB"]) only depends 
on the co-ordinates (£,r),J). A system of generalized action-angle co-ordinates in the 
sense of Nekhoroshev \2b^ is obtained by action-angle co-ordinates (0i, . . . , <pg, I\, . . . , Ig) 
of the reduced system and by (ip , . . . , ipg, J , . . . , J/) together with any sets of ^ 2m CT — 4 
symplectic co-ordinates on the Grafimannians G m<y 2 with m a > 3. 



Proof. The proof of this is similar to that of Lemma 13.61 The difference is that here we 
work in the Poisson structure instead of in the symplectic structure. While this saves us 
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the work of specifying a local symplectic structure on T*S e , the number |m -(m CT — 1) of 
global co-ordinates L ik exceeds the dimension 2m a — 4 of G m(T ,2 and the syzygy 

i<keI tT 

for m a > 4 is accompanied by increasingly more syzygies. The fact that J a has periodic 
flow with respect to the Dirac bracket can be proved as in Theorem 14.11 □ 

By the separation of variables in the next section we give a local symplectic co-ordinate 
system on T*S and are then able to compute I\, . . . , Ie explicitly. 

Theorem 4.7. (Energy-Casimir mapping, regular part) The fibre over each weakly reg- 
ular point of the Energy-Casimir mapping 

EC : T*R n+1 — > R i+2 

(x, y) ' y (H, J , ... , Ji) 

is given by 

T h>j = £C- 1 (h,j) = H- 1 (h)n |J J-'ifi). 

The structure of J-'hj is that of a double bundle in view of the two properties 

e 

(7=0 

where R is the reduction mapping of Theorem \4-5\ The fibre of each of the two mappings 
is a T £+1 bundle over the image of the other one, yielding the diagram 




Relative equilibria are given by the critical points of the amended potential V^: £ CT = for 
m a = 1 and £^ = j a /cu a , uj a = \Jb a — f3 otherwise, where (3 is determined from = 1. 
The corresponding critical value of the Energy-Casimir mapping is 

There are no other critical values of SC in the range 

{(h,j)eR i+1 :j a >0^ =o ^j} 
of regular reduction. 
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Proof. The double bundle structure appears because we chose to treat the actions J that 
come from the group G apart from the actions I that appear in the reduced system. In 
the end we obtain a T^ +1 -bundle, but for now we only understand the actions J because 
they are related to the Momentum mapping of G. In particular for the relative equilibria 
(relative to G\) this splitting is natural. 

Relative equilibria are solutions of VH = |/3VC*i + 5VC*2. Since C2 = the last n 
components of these equations give r\ = and 6 = 0. Hence V^H = (3£ needs to be 
satisfied. This leads to £ a b a — j^/^ — /3£<r = 0, which gives the result. The multiplier (3 
is determined by the constraint C\ = 1. □ 

Lemma 4.8. The critical energy h of the relative equilibria is a convex function of the 
momenta j a . 

Proof. The function h(j) is implicitly defined by h = f(j,(3) and g(h,j,(3) = 0, where / 
is given by the Energy-Casimir mapping and g is given by the constraint. We are going 
to prove that the Hessian of h(j) is positive semi-definite. By implicit differentiation we 
find 

d ja h = djj - dfsf^^- = u a + bv/ujv - f3/uj a = 2^/b a - (3 = 2u a 

where the identity dpf /dpg = (3 has been used. To compute the second derivative we 
therefore need the derivative of (3. It is given by 

d . 8 = -^ = -- 
jJ d fi g ^ ^ W 3 

Combining the two first derivatives gives the entries of the Hessian as 
8 2 h 2 
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This is a rank one matrix with one positive eigenvalue, and hence h(j) is a convex function. 

□ 

In order to describe the invariant manifolds of the invariant sets contained in lower dimen- 
sional invariant Neumann systems defined by j a = for some (or all) o singular reduction 
is used. 

Theorem 4.9. (Singular Reduction) The Hilbert mapping \ of the G-action is the direct 
product of the Hilbert mappings of the individual factors 0(m a ), as given in Lemma [X7 . 
The reduction mapping x '■ T*M n+1 — > M 3 ^ +1 ) is given by 

((^x; V a x), (V a y; V a y), (P CT x; V a y)) = (2V a , 2T a , S a ) 
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with syzygies 2V a T a = ~S„ whenever m a = 1. The reduced brackets are 

2V 2V 
{V a ,T T } = S T (5 aT - — {V a , S T } = 2V tT (5 rTT - —I), 
C-l Oi 

2V 2T S — 2T S 
{To, S T } = -2T a (S aT - — -), {T a ,T T } = — — , 

{S CT ,S T } = 0, {K,K} = 0. 

This bracket has rank 2£ and the £ + 3 Casimirs 

d = 2^K, £2 = $^ and ^ = 4^-5^. 

From the latter we obtain at nonzero values the £ + 1 actions J a = y/W a . The reduced 
Hamiltonian of the degenerate Neumann system reads 

e 

H(V,T,S) = J2T„ + b a V a . (18) 

(T=0 

Proof. The brackets with equal indices on the left hand side are found by direct calculation 
using (J3J). Only the additional terms need to be computed, for the [..,..] brackets see 
Lemma 13.71 Since C% = 2V and C2 = S the additional terms are also already known. 
The other brackets follow similarly by direct computation. □ 

Remark 4.10. For the co-ordinates with m a = 1 we may as well keep the original x a , y a 
and avoid the £ — £ syzygies. Note that the rank of the bracket remains 2£ with £ + 3 
Casimirs. All S a can be eliminated using the syzygies and fixing the values w a of the 
Casimirs W a , albeit at the cost of extracting a square root. The two sheets of ±S a are 
glued together smoothly at S a = 0. In addition using Ci, C 2 this projection of the reduced 
space is of dimension 2£, the same as the dimension of the regularly reduced phase space. 



5 The Rosochatius system 

Regular reduction of the degenerate Neumann system leads to an open subset of T*S l as 
described in Theorem 14.51 The reduced phase spaces J -1 (/i)/G M can be interpreted as the 
regular symplectic leaves of T*13L n+1 /G, which we embedded in R 3 ^+ 3 ) in Theorem 14.91 
The remaining symplectic leaves of T*M. n+1 /G can be given an interpretation of "filling 
up the remaining part" of T*S e , as already indicated in Remark 14.101 

Theorem 14.51 suggests to consider the following "unfolding" of the reduced Neumann 
system on T*S e C T*W e+1 with co-ordinates (£, rj) satisfying the Dirac brackets (T4|) in 
dimension 2£ + 2. The Hamiltonian is given by 

mv) = \i2vi+b„e + ^, (19) 

(7=0 
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where w a , o = 0, ...,£ are the unfolding parameters. When reducing the degenerate 
Neumann system we find w a = j% > 0, and it is this same condition that is used in [29| 
2H EE] to define the Rosochatius system on T*S l . Since b Q < . . . < b £ the symmetry group 
Z e + l = 0(1) x ... x 0(1) is discrete. 

Theorem 5.1. Singular reduction of the Z^ 1 -symmetry of the Rosochatius system leads 
to the same reduced dynamics as the singularly reduced degenerate Neumann system. The 
Hilbert mapping x '■ T*M. l+1 — > M 3 ^ +1 ) is given by the basic invariants 

2 ' u 2 2& 

which satisfy the Poisson bracket relations of Theorem \4-9\ and turn the Hamiltonian (j!8p 
into (USD- 

Proof. This follows again by direct computation, e.g. 

{Xt, Jt} = tyr^r ~ h r} a -pr{5 aT - — — J - — r/ r (d (JT - 



K = ^, T (T = ^ + ^, and S a = ^ Va 



2T CT • (S,- 2T t • S a 



□ 



In this way the Rosochatius system is a 2^-fold covering of the (singularly) reduced de- 
generate Neumann system. 

Remark 5.2. For coefficients b a of the degenerate Neumann system with m a = 1 we may 
keep the original x a ,y a variables ; for these indices w a = and the Rosochatius system 
is a 2^ +1 -fold covering of the resulting reduced system with G reduced instead of G. For 
coefficients b a of the degenerate Neumann system with m a = 2 we may similarly choose 
to reduce with respect to the 50(2)-factor instead of the 0(2)-factor. In particular, 
if m = . . . = mi = 2 then singular reduction of the degenerate Neumann system with 
respect to the subgroup £0(2) x. . .xSO(2) of G = 0(2) x . . . xO(2) yields the (unreduced) 
Rosochatius system, as had already been remarked in 



5.1 Separation 

In this section we consider the reduced system in its own right and separate the Hamil- 
tonian of the Rosochatius system in elliptical-spherical co-ordiates. Note that for the 
integrability the condition w a > is not needed. However, the dynamics with w a < is 
completely different as then there are solutions that blow up in finite time. 

Theorem 5.3. The Neumann system on S with additional potential \ ^2,w CT / ^ can be 
separated in elliptical- spherical co-ordinates Ui. The general solution of the Hamilton- 
Jacobi equations is 

S = - [ — -7— dz 
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where A(z) = Yl( z ~ ^o) an d the integral is defined on the real part of the hyperelliptic 
curve 

r = {(z, C) e C 2 : c 2 = -Q(p; z)A(z) + Q(b, w, z)} 

of genus i. The separation constants are pi in 
Q(p; z) = z E + 2p 1 z e ~ 1 + --- + 2p t 

where p\ is the value of the energy. The polynomial Q is given by 

i t 

Q(b,w; z) = ^Wa- Y\(z - b T ){b a - b T ) . 

cr=0 Tj^cr 

Proof. The separation proceeds as in the case without additional potential, see [23]. Start- 
ing from the Lagrangian in spherical co-ordinates defined by 

<T=0 

the co-ordinates are the roots Ui, i — 1, . . . , £ of the rational function f(z) = 0. Because 
the b a are different, and f(z) has a pole of first order at each b a the graph of f(z) shows 
that the roots u a satisfy 

bo < Ui < bi < ■ ■ ■ < ue < bi . 

The denominator of the rational function f(z) is the polynomial A(z). Define the poly- 
nomial in the numerator of f(z) as U(z) = Y\i z ~ u i)- Then we have 



e a = R**f(z)dz 



z=K " A'(b a ) ' 

From the standard approach, see [8J for the details, we know that after introducing the 
momentum pi conjugate to Ui the kinetic and original Neumann-potential energy read 



2 ^ u 2 

<T=0 1=1 

The additional potential V w in terms of u becomes 

A\K) 



\ '„■{ ») = ~ £ """ 

(7=0 
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so that 

H(u,p) = T(u,p) + V(u) + V w (u) (20) 

and the Hamilton-Jacobi equation reads H(u,dS/du) = p\. The general solution for S 
depends on I integration constants. We write the Hamilton-Jacobi equation in the form 
= ^2h(ui,pi)/U'(ui) where h(ui,pi) only depends on in and pi. The kinetic energy 
T(u,p) already has this form. To achieve this form for V{u) and pi the following two 
identities due to Jacobi are used, see Appendix [S] For any polynomial U of degree £ with 
roots Ui we have 



5> = E7tS 



i=i i=i 



If in addition we have a polynomial P(z) = p\z e 1 + ■ ■ ■ + pe with arbitrary coefficients 
then 



i=i y 1 



It remains to write V w {u) in the desired form. Consider the partial fraction decomposition 
of £~ 2 with respect to b a , 

Now exchange the summation over i in this formula and a in V w to find 

V 1 1 ^'(U 1 Q(b,w;ui) 

2^U'(u i )2^ W ' 7 b a -Ui j^U'{ Ui ) 2A( Ul ) ' ^ 

where 

i 

Q(b, w; m) = ]T — A(ui)A'(b a ) . 

O a Ui 

(7=0 

Hence, every term in (12"U|) is proportional to l/U'(ui) (up to the constant c = | an d 
if reads 

The Hamilton-Jacobi equation H(u,dS/du) = p\ becomes 



£ 1 

= Yl TpT~\ h ( Uh ds / du i) 



1 U'(u 

i=i v 
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where after shifting p\ by c we find 



h( Ui , Pi ) = -2A{ui)pt - X -u\ - P( Ui ) + 




To complete the proof define 

Q(p; z) = z e + 2P(z) =z e + 2p l z e ~ 1 + --- + 2p e , 



and the general solution of the Hamilton- Jacobi equation depending on £ separation 
constants p^, i = 1, . . . , £ is 



where ( is found by solving h(z, (/(2A(z))) = 0. To make the integral well-defined in the 
presence of the square root it is defined on the hyperelliptic curve given by 



Remark 5.4. The genus of the curve T is £, which is unchanged by introduction of w a 
since QA has degree 2£ + 1 while Q only has degree £. Hence for the reduced system the 
genus equals £, which is the number of different coefficients minus one. 

When £ = the curve is algebraic, this corresponds to the Neumann system with all 
axes equal, hence it is the geodesic flow on the sphere with all orbits closed, see Section FPl 

When £ = 1 the curve is elliptic and explicity given by (recall from reduction that w a 
is the value of J 2 ) 



If Jo = J\ = the Neumann system on the circle is recovered. 

When £ = 2 the curve is hyperelliptic, containing the classical Neumann system on 
the sphere special case for J = J\ = J 2 = 0. 

Remark 5.5. For £ = 1 the curve given in the previous remark is the same as for 
the Lagrange top. This equivalence is classical and has first been described by Klein Sz 
Sommerfeld pR)j : when the spherical Lagrange top is described on the double covering 
SU(2) = S 3 of SO (3) then the equations of motion of the Lagrange top are those of 
the Neumann system with potential ^(^o + x i) + ^i( x 2 + x V)-> hence £ — 1, mo = mi = 
2. As a result we can conclude that this degenerate Neumann system has Hamiltonian 
monodromy, inherited from the monodromy of the Lagrange top jl]. 

The relation between the Lagrange top on 5*0(3) and the degenerate Neumann system 
on S 3 is interesting because the action of the symmetry group T 2 in the Lagrange top is 
not free, although the constituting S'0(2)-actions are both free, while for the Neumann 
system each SO (2)-action already has the trivial fixed points. The resulting period lattice 
for the Lagrange top is generated by (Jo ± Ji)/2, and not by Jo, J\. This is an expression 
of the fact that SO(3) has half the volume of 5* 3 , which covers it twice. 




c 2 



Q(p; z)A(z) +Q(b,w;z) . 



□ 



C 2 = -(z + 2 Pl ){z - b )(z - h) + 4{z - 



b 1 ) + J 2 1 (z-b )(b 1 -b ). 
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5.2 Actions and Frequencies 

The degenerate Neumann system has up to i + 1 global actions J a , and for simplicity 
during the rest of this section we assume that this maximal number is realised, i.e. 1 = 1 
and all m a > 2. In addition to the I + 1 actions J a the separation of variables of the 
reduced system gives £ non-trivial actions ij. The integration proceeds between the branch 
points of the hyperelliptic curve T. Since the co-ordinates u cover only one 2 f+1 -tant of 
the sphere S the signs of the variables £, r\ need to be recovered. This is very similar 
to the non-degenerate Neumann system, see [8j for the details. Recovering the sign is 
crucial when some or all J a vanish because then the co-ordinate planes can be crossed 
and hence a passage to another 2^ +1 -tant takes place. The co-ordinate £ CT becomes zero 
when Ui = b n for some % G {1, . . . , £}. Consequently, when one of the endpoints of the 
segment [z2i-i, is equal to b a then the co-ordinate £ CT changes sign along the cycle q 
encircling this segment of positive ( in the complex z-plane. This means that the cycle ji 
on the Liouville torus corresponds to the cycle c$ taken twice. This proves the following 
result. 

Theorem 5.6. The non-trivial actions of the reduced degenerate Neumann system are 
given by 

where integration path 7, equals Ci if the segment [z2i-i, Z2i] contains none of bo, ... , be and 
7i = 2q otherwise. □ 

When the Neumann system is reduced by the discrete symmetry group of reflection of 
each co-ordinate then for the reduced system the non-trivial action is always given by the 
integral over q. 

Remark 5.7. The actions Ii are hyperelliptic integrals of the third kind, with poles at 
z = b a (recall that we assumed m G > 2 for all a) . In the non-degenerate Neumann system 
the actions are given by hyperelliptic integrals of the second kind [8]. Consider as the 
integration path a circle c a around the pole z = b a . Then 

1 / C 

Ja — <t , , \ d.Z . 

In this sense the hyperelliptic curve "knows" not only about but also about the trivial 
actions J a . The reason for this is the identity 

C 2 (M = Q(b, w; b a ) = w a l[(b a - b T f 

so that the residue of the action integrand &t z = b a is the action y/—w a . This is a 
common phenomenon: if there are simple actions (e.g. actions that are generators of a 
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global symmetry) then the non-trivial actions of the system are abelian integrals of the 
third kind and the integration around the poles gives the trivial actions or multiples 
thereof. 

The integrals Ji, . . . , 1% depend on the separation constants pi and on the trivial actions J G . 
Next to the values j = (jo, • • • , je) of these we also write h = (hi, . . . , hi) where hi is the 
value of the energy and hi are the values of the other smooth constants of motion Hi not 
related to the symmetry group G. These may e.g. be the values of i — 1 of the i + 1 
integrals ((?]), or the separation constants p 2 , • • • , Pi- The action mapping A then takes a 
regular value (h, j) eR f x W e+1 and assigns the actions A = (I, J). The period lattice and 
the frequencies are determined by the derivatives of A. The vector fields of the actions 
(written as a vector) can be expressed in terms of the vector fields of the 2£ + 1 functions 
(H, J) as 

Wlth T <tl 

The matrix T describes the period lattice of the integrable system. The frequencies of the 
Hamiltonian flows of Xjj, and Xj are denoted by O i? f2^ +(T+1 and taken as row vectors 
they form the frequency matrix Q. By definition the frequencies of a constant of motion 
are the coefficients of the decomposition of its vector field in terms of the vector fields of 
the actions, hence 



and therefore Q = T _1 and in particular 



so that the frequencies Qi of the Hamiltonian are the entries in the first row of Q = T _1 , 
hence Qi = T~*ei. 

In the present case the matrix T has a simple block structure caused by the trivial 
actions J a , namely 
/ di ar 

rp _ I OH dJ 

V id m 

Therefore, the frequency matrix is given by 
vaM- 1 f di n- 1 oi ^ 

V OH ) \dH> 8J 



n 



id 



€+1 



The (£+1) -dimensional identity matrix expresses the fact that J consists of global smooth 
actions, and hence each of their flows has frequency one. 

A particularly simple case occurs when the reduced system has one degree of freedom, 
I — 1. Then the upper left block is a scalar which is also the first entry in the frequency 
vector Qi. Frequency ratios are formed by dividing Qi by one of its entries, so the scalar 
dH/ di cancels and using the first entry the two frequency ratios are simply and J^- . 
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6 Global structure of the flow 



Let us first comment on some special properties that appear when the multiplicity m a of 
roots of the quadratic potential in the Neumann system is small. 

m a = 1. In this case the eigenvalue a u = b a is not equal to any of the other eigenvalues 
of the quadratic potential. The factor 0(1) of the symmetry group G is discrete and 
does not yield a component of the Momentum mapping J — we formally set J a = 0. 
Furthermore it may be preferable not to reduce with respect to this discrete factor, 
as this only leads to singular points without the reward of lowering the number of 
degrees of freedom. See also Remarks 14.101 and 15.21 

m a = 2. In this case reduction lowers the number of degrees of freedom by one, but 2 
equal eigenvalues of the potential do not yet contribute to a possible superintegra- 
bility of the degenerate Neumann system. In particular, this means that we may 
observe monodromy by only reducing the 5'0(2)-factor and not the whole 0(2)- 
factor. This keeps the sign in the value of the Casimir which we rather choose to 
be J CT = than J a = \/Lf k . In this way it becomes possible to find a closed 
non-contractible loop around a critical value of the Energy-Momentum mapping, 
which is located on the boundary of the image of the Energy-Casimir mapping. See 
also Remarks 13.51 14.44 EZl and 16.41 below. 

m a = 3. As soon as we have at least one group of three (or more) equal eigenvalues a u 
of the potential, the degenerate Neumann system becomes superintegrable. Still, 
as long as there are no groups of four or more equal eigenvalues the weakly regular 
values of the Momentum mapping J are in fact regular values — there are no syzygies 
constraining the L^. 

m a = 4. This case is still distinguished from the general m a > 5 in that the Grafimannian 
G^2 Q X (l I /At 1 1) Q 50(4) is a hypersurface of J~ 1 (| I), defined by the zero set 
of the second global Casimir Li 2 L 34: + L U L 23 + Li 3 L 2i (next to W a = L\ 2 + Lf 3 + 
L 2 A + L 2 23 + L 2 24 + L 2 3A ) of so (4). 

6.1 Integral— Casimir Mapping 

In order to describe another aspect of the reduced degenerate Neumann system we in- 
troduce yet another mapping JC in addition to EM. and EC considered in Section |H The 
name Integral-Casimir mapping captures the fact that there are two different kinds of 
conserved quantities in the system; those that come from the G-action and those that do 
not. The momenta of the G-action yield Casimirs for m a > 2. The other integrals (and of 
course the energy) are simply conserved quantities that are not related to the symmetry 
group G. 

Regular reduction together with separation of variables describes the system when 
all J a > and when the motion stays outside the singularities of the separating co- 
ordinate system. Singularities of the separating co-ordinate system occur when one of the 
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roots Ui is equal to bj, or when Ui = Ui + \ for some i. Hence, whenever the hyperelliptic 
curve has no double roots the motion on the corresponding torus is well described by 
the elliptical-spherical co-ordinates. To find the candidates where the Integral-Casimir 
mapping 2C does not have full rank we can therefore simply study the discriminant locus 
of the hyperelliptic curve. 

Lemma 6.1. The discrimiant locus for the reduced Neumann system with all w a 7^ is 
topologically a 2 e -tant without any internal structure. 

Proof. The position of the roots of the hyperelliptic curve is such that no root can be 
located at b a because setting z = b a in Q gives 



which implies that all the roots z = b a for w a = are not roots any more when w a ^ 0. 
For a real motion there must be a positive interval for each Ui between and 6j. 
Hence there must be two roots in each such interval. This uses up 2£ roots already. The 
remaining root must be to the left of bo, because ( 2 — > +oo for z — > — oo. Now all pairs 
of roots bounding the Ui intervals are separated by some b a , and hence for w a ^ the 
only possible collision is within a pair, leading to = Ui + \. Making k such pairs collide 
gives a corank k degenerate value in the image of the Energy-Momentum mapping. The 
corank k critical values give the (£ — /c) -dimensional faces of the 2^-tant. Making all pairs 
collide gives the corank i elliptic (relative) equilibrium point of the reduced system. This 
corresponds to the origin of the 2 £ -tant. □ 

For I = 2 the discriminant locus of the curve of the reduced system is explicitly given by 



When w a = then in addition the straight line given by Q(b a ) = needs to be added to 
the discriminant locus. 

6.2 Convexity 

The simple discrimiant loci described in the previous Lemma 16.11 are the fibres over J— 
space. In this view the J-space is the parameter space for the reduced system. Instead 
of the discriminant locus of the reduced system we now study the J-space as the image 
of the Casimir mapping of the full system. More precisely, we consider the image of the 
Momentum mapping of the T^ +1 -action on phase space given by J, for short the Casimir 
mapping. To obtain a compact image we fix the energy, and study the set of possible 
values j of J. 



Q{K) 
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For each j there is a unique elliptic (relative) equilibrium of the reduced system. This 
corresponds to the origin of the 2 £ -tant of Lemma [67TI and there is a corresponding energy. 
Fixing instead the energy and asking for all equilibria with any value j gives the boundary 
of the image of the Casimir mapping restricted to constant energy. 

Theorem 6.2. The image of the Energy- Casimir mapping is a convex set for sufficiently 
large h. 

Proof. It is important to use the actions J a as co-ordinates in the image, not their squares. 
The result follows from the explicit parametrisation of the discriminant locus. We equate 
— QA + Q to — {z — s) 2 (z — ri) . . . (z — r 2 e-\) and solve for j a and p^. This leads to 

^ i ^ n-\ 

(7=0 fc=l 

and 



Ja= A'(b a ) 




,fe=i 



Consider now a parametrisation of the corank £ stratum of the discriminant locus S, 
which gives the relative equilibria which mark the boundary of the image of the Casimir 
mapping. Hence we parametrise by the t double roots, s&, h — 1, . . . , I and find as a special 
case of the above formula the values where are made pairwise equal and r denotes the 
remaining single root 



and at the same time p = p(sk,r). Since p\ is linear in r we can eliminate r in favour 
of the value h — p\ of the Hamiltonian and the sum of the Sk- Also we introduce the 
symmetric functions tk of the Sk and find 

]* = -^^(bi + b e -H 1 + --- + t f X with u 2 a = h + b a -2t 1 . 

The parametrisation is linear in tk for all k > 1. 

For sufficiently large energy the boundary of the image of the Energy-Casimir mapping 
is given by the critical values j a (h,ti, . . . ,tt). This implicitly defines h c (j , . . . ,je), the 
value of the critical energy as a function of the momenta. By implicit differentiation we 
compute the gradient of the energy 

dh c 
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and the Hessian as 

J (7 J T I U a U T t=() 

This is a rank 1 symmetric matrix with nonzero eigenvalue 

whose eigenvector is (1/coq, . . . , I /cog). As a result the Taylor expansion of the critical 
energy is 

£ / £ \ ^ 

Mj + A) = 2j> CT A CT + 2^(^^j +0(A 3 ). 

If P > then the Hessian is positive semi-definite, and the boundary of the image of 
the Energy-Casimir mapping is the graph of a convex function, and hence the image is a 
convex set. Now 6j < s i+1 < b i+ ±, and with t\ = — s i we find 

i 

B -b < -ti < B -b e , B = ^b a 

cr=0 

Since all u„ must be real, we conclude that the image is convex if the energy is sufficiently 
large, namely when h > —2(5 — be) — bo. Rewriting ti in the roots gives 

t i 

V = Y[(h + Si + 2S), S = J2 S *- 
i=i i=i 

From the above /1 + 2S" > -2(B-b e ) -b + 2B-2b = 2b e -3b and therefore h + s k + 2S > 
2be - 2b > and P > as desired. □ 

Theorem 6.3. For sufficiently large fixed energy the image of the Casimir mapping J 
given by the T e+l -action is convex. Fork — > oo the boundary tends to a convex polyhedron. 

Proof. Since we have shown that h c (j , . . . , j a ) is a convex function, and h c is the boundary 
of the image of the Energy-Casimir mapping for h > —2(B—be)—b , every slice of constant 
energy through this convex set is again a convex set. This set is the image of the Casimir 
mapping applied to the surface of constant energy. In the limit h — > oo it becomes a 
convex polyhedron when rescaled by \fh. This amounts to setting u a = 1 in the limit, 
and so the non-linearity disappears. Even for finite h the boundary has zero curvature 
for I > 2. In fact the boundary is a ruled surface for £ = 2 and a higher dimensional 
analogue with an £ — 1 dimensional plane attached at each point of the hypersurface. □ 
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Remark 6.4. The above statement about convexity becomes particularly nice in the case 
that all m a = 2. Then the J a can all be chosen to be signed quantities, and are defined 
also for J a = 0, so that there is a global T £+1 -action on phase space. In this way the 
convex set extends through the faces J a = 0. Note that this convexity result is obtained 
even though the phase space is not compact. The compactness is achieved by restricting 
to the compact energy surfaces h = const. 
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A Jacobi Trick 

Let 

e 

U(z) = Y[(z- Ul ) 
i=i 

be a polynomial of degree £ with distinct roots «j. Then the derivative of U evaluated at 
the root Uk is 

e 

U'{u k ) = Yl ( u k - u i) ■ 

This gives a way to exclude one factor in U and thus the partial fraction decomposition 
of l/U is 

1 1 1 

— = Y - . 

U(z) ttU'( Ui )(z- Ui ) 
The identity (the use of which is sometimes called Jacobi's trick) 



^ Vim) 
i=i i=i v l > 

follows from calculation of 

z l 

• dz 



2tt J U(z) 
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along a countour that is a sufficiently large circle that encloses all poles Ui. The left hand 
side is found by calculation of the residue at infinity, while the right hands side is found 
from the sum over the finite residues. In a similar way the identity 

i = V-^— 

can be obtained. Lowering the exponent further the residue at infinity vanishes. Jacobi's 
trick is used in the separation of variables in order to introduce the separation constants. 
For powers of z greater than £ higher order symmetric functions of the roots are found. 
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